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A Coulomb lattice-gas model with a host lattice screening mechanism is adapted to describe the ordering 
phenomena in alkali-metal fullerides of body-centered-cubic structure. It is assumed that the electric charge 
of an alkali ion residing at a tetrahedral interstitial site is completely screened by its first-neighbor Ceo 
molecules. The electronic energy of the C^g ion is also taken into consideration as a charged spherical shell. 
By means of these assumptions an effective (short range) pair interaction is obtained between two alkali ions. 
The resultant lattice-gas model is analyzed by using two- and six-sublattice mean-field approximations. The 
thermodynamic properties are summarized in phase diagrams for different shell radii. 

PACS numbers: 64.60. Cn, 05.70.-a, 61.43. Bn 
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I. INTRODUCTION 

In a previous paper |Q] (henceforth referred to as I) 
we have introduced a lattice-gas model to describe the 
ordering of alkahne atoms intercalated into solid Cgo of 
face-centered-cubic (FCC) structure. This model now 
is adapted to a body-centered-cubic (BCC) host lattice 
because these structures are also observed in experiments 
(see reviews in refs. 

Lattice-gas (Ising) models provide an adequate de- 
scription of ordering phenomena in intercalated alloys. 
Difficulties arise, however, when determining the effec- 
tive pair interaction between two intercalated particles. 
In A^Ceo alkali- fullerides the alkali atoms transfer their 
s electron to Cgo molecules. As a result the Coulomb en- 
ergy becomes dominant in the formation of different or- 
dered structures. |^,|| At the same time these transferred 
electrons take part in the screening of Coulomb interac- 
tion between two alkaline ions. In the knowledge of the 
screening mechanism one could derive an effective pair 
interaction expected to be a short range one. Unfortu- 
nately, we have no satisfactory approach for determining 
the screening in these materials. 

In the present model we assume a simple screening 
mechanism. First of all the effect of dielectric media 
is considered viaa dielectric constant s. Furthermore, 
the charge of each alkali ion is completely screened by 
the first neighbor Cgo molecules. This means a uniform 
distribution of transferred s electrons to the first neigh- 
bor Cgo molecules. This simplification seems to be also 
an adequate approximation in a metallic system if the 
Thomas-Fermi length is not greater than the shortest A- 
Ceo distance (« 6.3 A ). 

In accordance with the above assumption the charge 
assembled on a given Cgo molecule may not always the 
same: its value depends on the number of alkali ions 



residing around the molecule. A series of UHF-MNDO 
calculations has confirmed that the electronic energy of 
a solitary Cqq ion may be approximated as E{x) — 



- bx - 



where the quadratic term mediates an inter- 



action between those alkali ions transferring charges to 
the same Cgo molecule(s). This contribution may be con- 
sidered as the electrostatic energy of a charged spherical 
shell with a radius R = 4.45 A because the term a + bx 
does not affect the thermodynamic behavior. The effect 
of dielectric media is taken into account with the appli- 
cation of the Kirkwood-Onsager theory as in I. Rabe 
et al. suggested such a spherical shell model choosing 
i? = 3.5 A (and £ = 1) which is equivalent to the shell 
radius formed by the carbon nucleus of Ceo molecule. ||] 
We use the concept of this spherical shell model because 
of its descriptiveness, i.e. the electronic energy of the 
Cg(^ ion will be estimated as E{x) — {xe)^ /2eR. 

Without repeating the formalism introduced in I, we 
now concentrate on the mean-field analysis of the model 
using six- and two-sublattice approximations. The re- 
sults are summarized in phase diagrams calculated for 
three different values of R. In the first case R is esti- 
mated by using the data of UHF-MNDO calculations. 
Accepting the suggestion of Rabe et al. , in the second 
case, we have determined the phase diagram for i? = 3.5 
A. Finally we have chosen an intermediate value for R to 
display a phase diagram different from previous ones. 

In the analysis of the above model we restrict our- 
selves to the rigid BCC structure. Our main purpose is 
to analyze the general features of the present model in 
the parameter region relevant for alkali-metal fullerides. 
A more rigorous comparison of the states appearing in 
FCC, BCC and body-centered-tetragonal (BCT) host 
lattices [^-10| goes beyond the scope of the present work. 
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II. THE MODEL 



In BCC structures there are only tetrahedral intersti- 
tial sites for the intercalated alkali ions. To describe the 
distribution of A"*" ions we introduce a site variable rji 
which is 1 if the interstitial site i is occupied by an alkali 
ion and for empty sites. The energy for any configura- 
tion of alkali ions is given by the lattice-gas Hamiltonian 



(1) 



characterized by the site energy e and the effective pair 
interaction Vij. Both quantities are determined by as- 
suming the screening mechanism mentioned above. 

The calculation detailed in I results in an effective 
pair interaction dependent on the ion-ion distance rtj as 
shown in Fig. |^ for R = 4.45 A. Henceforth the energy is 
measured in units of e^/(ea) where a is the lattice con- 
stant. Notice that Vij may have two values for certain 
values of r^. This effective interaction is repulsive for 
short distances and becomes attractive when increasing 
the ion-ion distance. In agreement with expectations the 
strength of Vij decreases rapidly with Vij . 




FIG. 1. Effective pair interaction vs. ion-ion distance for 
R = 4.45 A. 

The interaction mediated by Cqq ion(s) gives a repul- 
sive contribution to the total pair interaction until rij 
exceeds a threshold value. Obviously, the strength of 
this contribution increases when decreasing R. 

In the present model the site energy may easily be 
determined: 



-1.3752- + — 

ea seR 



(2) 



In the knowledge of the model parameters we are able 
to determine the energy of each ordered particle distri- 
bution. For this purpose the tetrahedral interstitial sites 
are divided into 12 simple cubic sublattices indicated by 
hexadecimal codes from 1 to C as illustrated in Fig. 2. 




FIG. 2. Tetrahedral interstital sites of BCC lattice are di- 
vided into 12 sublattices. 

Within this sublattice division we can distinguish 2^^ 
ordered states in which the sites belonging to the same 
sublattice are uniformly occupied or empty. Figure 3 
shows the energies per Ceo of these states as a function 
of alkali content (x). Here it is worth mentioning that the 
charges of Cgo molecules are uniform within an ordered 
state. Consequently, the energies plotted in Fig. 3 are 
the sums of the corresponding Madelung energy and the 
estimated electronic energy of Cgo- 




FIG. 3. Energy of 12-sublattice ordered states as a function 
of alkali- metal content for R = 4.45 A. Open diamonds and 
squares refer to states described by the two- and six-sublattice 
formalisms. The solid line connects the stable phases. 

Figure 3 demonstrates clearly that this model has only 
three stable ordered structures if i? = 4.45 A. For the 
stable a; = (6) state all the tetrahedral sites are empty 
(occupied). In the half-filled system {x — 3) the stable 
state is equivalent to the well known A15 structure, i.e. 
the sublattices 2, 4, 6, 7, 9 and B are completely occupied 
whereas the rest sites are empty. If the equivalent sub- 
lattices are combined the formation of this structure may 
be well described by using a much simpler two-sublattice 
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formalism. 

The states with the minimum energy for x = 1,2,4 
and 5 exhibit additional symmetry, more precisely, the 
occupations in sublattices a and a + 6 are equivalent for 
a = 1,...,6. The combination of the equivalent sublat- 
tices results in the six-sublattice formalism studied pre- 
viously. In this case the tetrahedral sites are divided 
into six interpenetrating BCC sublattices therefore the 
translation symmetry of the host lattice is conserved in 
the mentioned states. The role of these states becomes 
important when studying the system for lower R. 

The decrease of R gives an x-dependent contribution 
to the energies (see Fig. 3) without modifying the energy 
differences for fixed x. The minimum energy states for 
X — 1 and 5 become stable if i?/a < 0.3411. The simulta- 
neous appearance of the stable states for a; = 1 and 5 is a 
consequence of the particle-hole symmetry. In the x = 1 
stable state only the sites of one of the six sublattices 
are occupied. If R/a < 0.3342 then the system exhibits 
stable states for x = 2 and 4. In the x — 2 stable state 
the simultaneous occupation of the first neighbor sites is 
excluded. 

The above ground state investigations support that the 
subsequent mean-field analysis may be restricted to the 
two- and six-sublattice approximations. 

III. MEAN-FIELD APPROXIMATION 

The six-sublattice mean-field approximation was pre- 
viously introduced to describe the ordering processes in 
superionic Agl where the I~ ions form a BCC cage lat- 
tice for the mobile Ag+ ions. |0 This sublattice division 
is obvious because the number of tetrahedral interstitial 
sites is just six times greater than the points in the BCC 
host lattice. In this approach the equilibrium states are 
characterized by the average sublattice occupations cFi, 
{u = 1, ...,6). The energy per Cgo molecule is expressed 
by these quantities as 

7i = ^ £(7^ + i ^ J^r^iyCfr (3) 

where e is defined by Eq. and the mean-field coupling 
constants J^^t summarize the effective interactions {Vij) 
between a site of sublattice v and all of the sites belonging 
to sublattice r. Due to the symmetries of this sublattice 
structure jl^] there are only three independent mean-field 
coupling constants. The numerical calculations yield 

Jii - ... = Jm = -1.16445— , 

ea 

e2 

Ji2 = ... = J65 = +0.08 2 51— , (4) 

ea 

Ji4 = ... = J63 = -0.30 4 68— , 

ea 

if i? = 4.45 A. 



In the second approach we distinguish only two sub- 
lattices (a and /3) as described in the previous section. 
Now the mean-field energy per Cgo molecule is defined 
by the following expression: 

7i = S^ecr^ + - ^ J;^3CrrCrs (5) 
r r,s 

where r, s = a, (3 and the values of J^^ are determined 
numerically for the same radius: 

JL. = 4/3 = -1-83360^ , 

4/3 = J'pc = +0-69451^ . (6) 

Within this formalism the twofold degenerate A15 struc- 
ture is given as a a = {1 — a^) = or 1. Obviously, the 
energies of the empty (ctq — <jfj — Q) and completely 
occupied {a a = = 1) states are equivalent to those 
suggested by the six-sublattice formalism. 

Using the same method one can also evaluate the 
mean-field coupling constants of the 12-sublattice ap- 
proximation. The values of these parameters are, of 
course, strongly related to the previous ones defined by 
Eqs. and (^). More precisely, Jyr and J'^^ are the 
linear combinations of the mean-field coupling constants 
of the 12-sublattice approximation. 

The equilibrium state characterized by sublattice occu- 
pations is determined by minimizing the Gibbs potential 
for fixed temperature and chemical potential. Following 
I, the results of this mean-field analysis are summarized 
in phase diagrams. 

In agreement with the particle-hole symmetry the 
phase diagrams are also symmetric, as illustrated in 
Fig. 4. In this figure the dashed line represents the critical 
temperature of the continuous sublattice ordering from 
the high temperature {a a — o'fs) distribution towards the 
A15 structure. Below the critical temperature aa differs 
from and the alkali content per Cgo molecule is given 
as a; = 3((Tq + ct/?). In the subsequent phase diagrams 
this state is denoted as A referring to the symmetry of 
A15 structure. 
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FIG. 4. Phase diagram suggested by mean-field approxima- 
tion for R — 4.45 A . The dashed line represents a continuous 
two-sublattice ordering process. The phases are indicated by 
the alkali content of ordered states. 

The two-sublattice ordering process is analogous to 
the formation of antiferromagnetic order in Ising mod- 
els. The critical (Neel) temperature may be expressed 
by the coupling constants as a frmction of that is, 



ksTN = [J'ali - J'aa) 



(7) 



This transition temperature has a maximum at a; = 3 and 
vanishes when x goes to or 6. Accepting the numerical 
value e = 9 deduced by Sanguinetti and Benedek 
from the charge-induced vibrational shift the above ex- 
pression predicts Tjv « lOOOK for a: = 3. Instead of this 
continuous transition, at low temperatures the system 
segregates into two phases with different alkali contents 
(indicated by solid lines in Fig. 4). The values refer to 
the alkali contents in the ordered states at T = 0. 

We have determined the phase diagram for i? = 3.5 A 
{R/a — 0.303) suggested by Rabe et al. In this case 
only half (3 < a: < 6) of the symmetric phase diagram is 
represented in Fig. 5. 
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FIG. 5. Phase diagram in mean-field approximation for 
R — 3.5 A exhibits a eutectoid phase transition at a; ~ 4.5 . 

This figure shows stable ordered structures for a; — 3,4, 
5, and 6 at low temperatures. The a; = 3 state transforms 
continuously to the random distribution with a critical 
temperature given by Eq. (|^). A part of this curve is 
denoted by the dashed line in Fig. 5. As demonstrated, 
state A can exist for a wide range of x. The ordered 
structures with nominal compositions x — A and 5 be- 
come unstable when increasing the temperature above 
kBT4 = 0.1835eV(ea) and fc^Ts = 0.1971eV(ea). Fur- 
thermore the model exhibits a eutectoid phase transition 
at a; ~ 4.5. This phenomenon implies the stability of an 
ordered structure with x ~ 4.5 for smaller R. 

Here it is worth mentioning that the cubic symmetry 
is broken in the ordered structures if a; = 1, 2, 4 and 
5. As a consequence, the diffusion becomes anisotropic if 



the particles (holes) reside predominantly in two (second 
neighbor) sublattices. It seems reasonable to suggest 
that a distortion of the host lattice occurs because of the 
breaking of cubic symmetry related to the ordering pro- 
cess. A similar phenomenon is investigated by O'SuUivan 
et al. with MD simulations when studying the ordering 
process in the superionic phase of Agl. 

As mentioned above there exists a very narrow range 
of R where the ordered state a; = 5 is observable whereas 
the state x = 4 is unstable. In order to illustrate typical 
behavior within this region we have evaluated the phase 
diagram for R/a — 0.337. 

In Fig. 6 the dashed line refers to the critical tem- 
perature of the continuous sublattice ordering described 
above. In this phase diagram there is a region denoted 
as A-l-A' where two A phases coexist with different com- 
positions. At lower temperatures (see region A+6) the 
difference between CTq, and ap disappears for large alkali 
content (a: w 5.5). State 5 (with five of six occupied 
sublattices) becomes stable if ksT < O.I6I5e^/(£a). 
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FIG. 6. Phase diagram for R/a = 0.337 . 

The above series of phase diagrams visualizes the elec- 
tronic energy effect of the Cqq ion on the thermodynamic 
behavior. In contrast to the previous experiments and 
calculations the present model suggests the appearance 
of state 5 if state 4 is stable. Very recently, however, Lof 
et al. 15 1 have found a local minimum near x — 5 when 
measuring the temperature and concentration dependent 
conductivity of potassium doped Cgo films. This observa- 
tion may be interpreted as an experimental indication for 
the existence of such a state. Here it is emphasized, that 
this state is sixfold degenerated therefore a polydomain 
structure is expected to appear in experiments. 



IV. SUMMARY AND CONCLUSIONS 

Using two- and six-sublattice mean-field approxima- 
tions we have studied a lattice-gas model of intercalation 
alloys in which the Coulomb interaction between the in- 
tercalated elements is screened out by the charges dis- 
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tributed on the BCC host lattice. We have assumed that 
the intercalated particle residing in a tetrahedral inter- 
stitial position transfers e/4 charges to the (four) near- 
est neighbor host lattice molecules. The energy of the 
charged host lattice molecule is treated as the electro- 
static energy of a charged spherical shell of radius R. The 
family of alkali-fuUerides is the best candidate to satisfy 
the above conditions. This fact has motivated the choice 
of parameters when determining the thermodynamic be- 
havior suggested by this model. Although our analysis is 
restricted to alkali-fuUerides we can derive predictions for 
earth alkalines if 2e is sunstituted for e in all the formulae 
including the energy unit in the figures. 

The model takes the electrostatic interactions exactly 
into consideration at zero temperature, i.e. it reproduces 
the Madelung energies for the ordered structures. Owing 
to its simplicity the model has only two parameters: the 
lattice constant is taken from experiments, the value of 
R is estimated by a series of UHF-MNDO calculations. 
In this case the model suggests only three stable ordered 
structures. Besides the empty {x = 0) and fully occu- 
pied (x = 6) states we have found a stable A15 structure 
{x = 3). The fully occupied {x = 6) state is observed 
experimentally for A=K, Rb and Cs. For smaller alkali 
elements (Li, Na), however, the interstitial voids are too 
large compared with ionic sizes. Yildirim et at suggest 
clusters of four to nine sodium atoms inside the octa- 
hedral voids of the FCC lattice, jl^] The structure of 
NagCeo exhibits an FCC host lattice with single occu- 
pation of tetrahedral sites and fourfold filling of the oc- 
tahedral voids. Obviously, the description of these 
former structures goes beyond the validity of the present 
lattice-gas formalism. 

According to the electrostatic force calculations the 
FCC structure is preferable to the BCC one for x < 3. 
1 4, 5 1 This theoretical prediction agrees with experiments. 
|2,3| Up till now only the BasCgo compound has been 
found to exhibit the A15 structure. 

In the present model the robust A (A15) state can ex- 
ist in a wide range of temperature and alkali content as 
indicated in the phase diagrams. On the analogy of an- 
tiferromagnetic ordering this particle arrangement trans- 
forms continuously to the random distribution when the 
temperature is increased. 

The most surprizing prediction of the present calcula- 
tion is the existence of the A5C60 phase when the A4C60 
compound is stable. In these ordered structures the cu- 
bic symmetry is broken therefore lattice distortions are 
expected to accompany the ordering processes. For ex- 
ample, the A4C60 compounds have been observed with 
BCT structure wich may be considered as a distorted 
BCC lattice with four of six sublattices to be occupied. 
The present model permits the existence of the stable 
A4C60 compounds together with A5C60 ii R/a < 0.334. 
The X = 5 composition, however, has never been detected 
in experiments. The breaking of cubic symmetry in these 
former states raises many questions related to the role of 
lattice distortion. The investigation of the present model 



making use of a BCT structure is in progress. The results 
of this analysis are needed to enable a complete picture 
to be drawn. 
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